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Abstract. The integrated galactic initial mass function (IGIMF) is computed from 
the combination of the stellar initial mass function (IMF) and the embedded cluster 
mass function, described by a power law with index /3. The result of the combination is 
a time-varying IMF which depends on the star formation rate. We applied the IGIMF 
formalism to a chemical evolution model for the solar neighbourhood and compared 
the results obtained by assuming three possible values for /3 with the ones obtained 
by means of a standard, well-tested, constant IMF. In general, a lower absolute value 
of /3 implies a flatter IGIMF, hence a larger number of massive stars, higher Type la 
and II supernova rates, higher mass ejection rates and higher [a/Fe] values at a given 
metallicity. Our suggested fiducial value for /3 is 2, since with this value we can account 
for most of the local observables. We discuss our results in a broader perspective, with 
some implications regarding the possible universality of the IMF and the importance of 
the star formation threshold. 



1. Introduction 

The initial stellar mass function (IMF) is of primary importance in galactic chemical 
evolution models. The IMF regulates the relative fractions of stars of different masses, 
hence their relative contribution to the chemical enrichment of the interstellar medium 
(ISM) is tightly related to this quantity. For this reason, the analysis of abundance ratios 
in galaxies may allow one to put robust constraints on both the normalization and the 
slope of the IMF (Recchi et al. 2009; Calura et al. 2010). 

The Solar Neighbourhood (S. N. hereinafter) can be considered the most valuable envi- 
ronment to achieve constraints on the main parameters regulationg chemical evolution 
models, since it is definitely the best studied Galactic environment and many obser- 
vational investigations devoted to its study provide us with a large set of observables 
against which models can be tested. These observables include diagrams of abundance 
ratios versus metallicity, which are particularly useful when they involve two elements 
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synthesised by stars on different timescales. An example is the [or/Fe] vs [Fe/H] dia- 
gram, since a elements are produced mostly by massive stars {m > SMq) on very short 
(< 0.03 Gyr) timescales, while type la supemovae (SNe) produce most of the Fe on 
timescales ranging from 0.03 Gyr up to one Hubble time (Matteucci 2001). This diag- 
nostic is a strong function of the IMF, but depends also on the assumed star formation 
(SF) history (Matteucci 2001; Calura et al. 2009). Another fundamental constraint is 
the stellar metallicity distribution (SMD), which depends mainly on the IMF and on 
the infall history (hence on the star formation history) of the studied system. Another 
example of a useful diagnostic test for the IMF is the present-day mass function, i.e. the 
mass function of living stars observed now in the Solar Vicinity (Elmegreen & Scalo 
2006). 

The integrated galactic initial mass function (IGIMF) originates from the combination 
of the stellar IMF within each star cluster and of the embedded cluster mass function 
(CMF). It relies on the observational evidence that small clusters are more numerous in 
galaxies and that the most massive stars tend to form preferentially in massive clusters 
(Weidner & Kroupa 2006). The IGIMF is star-formation dependent, hence it is time- 
dependent and its evolution with time is sensitive to the star formation history of the 
environment. 

In this paper, we use all the local observables to study the IGIMF and its effects on the 
chemical evolution of the solar neighbourhood. The results obtained with the IGIMF 
are compared to those obtained with a non-star-formation dependent (hence constant in 
time), fiducial IMF. The aim is to derive some contraints on the main unknown param- 
eter of the IGIMF, i.e. the index y6 of the power law expressing the embedded CMF. 
This paper is organized as follows. In Section 2 we present a description of the theoret- 
ical scenario behind the IGIMF. In Section 3 we describe the chemical evolution model 
for the Solar Neighbourhood. In Sect. 4 we present our results and finally in Sect. 5 
some open problems regarding the IMF are discussed and some conclusions are drawn. 



Table 1 . Solar neighbourhood observables, parameters to which they are most 
sensititve and references. 



Observable 


Parameter 


Reference 


SFR Surface density 
type la SNR 
type n SNR 
Gas surface density 

Stellar surface density 
Stellar abundance ratios 
Stellar Metallicity distribution 
Present-day mass function 


SF efhcency 

Integrated SF history, IMF 
SF efficiency, IMF 
SF efficiency, IMF 

SF history, IMF 
SF history, IMF 
SF history, IMF 
SF history, IMF 


Rana (1991) 
Cappellaro (1996) 
Cappellaro (1996) 
Kulkarni & Heiles (1987) 
Olhng & Merrifield (2001) 
Weber & de Boer (2009) 
various authors 
Jorgensen (2000) 
Miller & Scalo (1979) 
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Figure 1. Left panel: IGIMFs for different cluster mass functions. Upper panel: 
yS = 1; central panel: fi -2; lower panel: p - 2.35. In each panel we have considered 
7 possible values of SFRs, ranging from lO"'* M^yr^^ (lowermost solid lines) to 
100 M^yr^^ (uppermost solid lines), equally spaced in logarithm. Right panel: the 
open circles are the "fitness" as a function of (3 for various models with different SF 
efficiencies (indicated by the numbers beside each open circle in units of Gyr '.). 
The horizontal line indicates the fitness value computed for the standard model. 

2. The (integrated galactic) initial mass function 

The main equation used to calculate the IGIMF is (see the contributions by Pflamm- 
Altenburg et al. and Weidner et al.): 

^lGlMF(m;iJ/(t)) = I ^C/W < '«max)^ecl(Mecl)<iMed, (1) 

where ifr is the star formation rate (SFR). The canonical stellar IMF is ^(m) = km~", 
with 0- = 1.3 for 0.1 Mo < m < 0.5 Mq and a ^ 2.35 for 0.5 Mq < m < mmax- The 
upper mass m^ax is a function of the mass of the embedded cluster Meci: this is logical 
if one considers that small clusters do not have enough mass to produce very massive 
stars. Star clusters are also apparently distributed according to a single-slope power 
law, ^eci °^ (Lada & Lada 2003). In this work we have assumed 3 possible values 
ofyS: l,2andZ35. 

Meci,min and Med.maxCiACO) are the minimum and maximum possible masses of the clus- 
ters in a population of clusters, respectively, and wJmax = 'WmaxC^eci)- For Meci,min we 
take 5 Mq (the mass of a Taurus-Auriga aggregate, which is arguably the smallest star- 
forming "cluster" known). The upper mass of the cluster population depends instead 
on the SFR and that makes the whole IGIMF dependent on ifj. 
The standard IMF is a two-slope power law, defined in number as: 

^ , , f 0.19 -m-^-^^ ifm<2Mo 
6.^(m) = I Q 24 . ifm>2Mo, 



This equation represents a simplified two-slope approximation of the actual Scalo 
(1986) IMF. The IMF and all the IGIMFs are normalised in mass to unity as the standard 
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Figure 2. Upper panels: SFR, type la and II SN rates, stellar and gas density 
versus time (left), abundance ratios (center) and stellar metallicity distribution (right) 
calculated with /? = 1. Middle panels: same as above but with [3-2. Lower panels: 
same as above but with fi - 2.35 (in the left panels only the case with SF efficiency 
1.5 Gyr"' is shown). For details see Calura et al. (2010). 



IMF: 

J" m^{m)dm - \. (3) 

In Fig. [U we show the IGIMF as a function of the SFR for the three values of 
P considered in this work, compared to our standard IMF. In general, a lower value 
of p implies a flatter IGIMF, and a hence higher relative fraction of stars with masses 
m > \Mq. 



3. The chemical evolution model for the Solar Neighbourhood 

The main chemical evolution model is described in detail in Calura et al. (2010). The 
model is calibrated in order to reproduce a large set of observational constraints for the 
Milky Way galaxy (Chiappini et al. 2001). The Galactic disc is approximated by several 
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independent rings, 2 kpc wide, without exchange of matter between them. The Milky 
Way is assumed to form as a result of two main infall episodes. During the first episode, 
the halo and the thick disc are formed. During the second episode, a slower infall of 
external gas forms the thin disc with the gas accumulating faster in the inner than in the 
outer region ("inside-out" scenario, Matteucci & Fran9ois 1989). The process of disc 
formation is much longer than the halo and bulge formation, with time scales varying 
from ~ 2 Gyr in the inner disc to ~ 7 Gyr in the solar region and up to 20 Gyr in the 
outermost disc. In this paper, we are interested in the effects of a time-variyng IMF in 
the Solar Neighbourhood. For this purpose, we focus on a ring located at 8 kpc from the 
Galactic centre, 2 kpc wide. The model includes the contributions of type la SNe, type 
II SNe and low and intermediate mass stars to the chemical enrichment of the ISM. A 
SF threshold is adopted (cr,/, ~ IMqPc'^) to reproduce various local features, including 
abundance gradients (Colavitti et al. 2009). The IGIMF is allowed to vary as a function 
of the SFR, which in turn is a function of cosmic time. The IGIMF is calculated as a 
function of the SFR according to Eq. [T] In Table [T] we show the solar neighbourhood 
observables used in this paper, with the main parameters on which they depend and 
references. 

In order to quantitatively compare our results with the observables considered in this 
paper, we define the fitness quantity as: 

fUness = ^;S = E,- >-(0[»bs(0 - theo(/)]^ 
\ + 6 max{[obs(/)]2, [theo(0]2) 

where, for the i-th value of each considered parameter, obs (/) and theo (/) are the 
observed values and the predictions of the model, respectively. The weight w{i) is used 
to give each set of observables the same statistical weight. 

The closer fitness is to 1 , the better the model is in reproducing the observations. In 
the right panel of Fig. [T] we show the "fitness" as a function of /? for all the models 
considered in this paper (see Sect. H]). 



4. Results 



Some results of our study are summarized in Fig. [2] Here we show the results of 
the comparison between observables and model predictions computed with different 
assumptions for the cluster mass function index yS. In general, the local observed quan- 
tities are plotted with their error bars. For a given value of /3, the curves of different 
types are the model results computed with different assumptions for the SF effiency, 
reported in the legends. From the time evolution of the SFR, type la and II supernova 
rate (SNR), gas and stellar mass density (left panels) computed with various assump- 
tions for the parameter /?, we can see that lower values of y6 imply higher SNRs, higher 
gas mass densities and in general lower mass locked up in living stars and remnants. 
This is visible from the comparison of the results computed with the same values for 
the SF efficiencies and different values for /?. Another important issue regards the de- 
pendence on the star formation threshold: some values of p (fi = I, fS = 2) show SF 
histories indendent from the SF threshold. This is basically due to the large mass re- 
turned by dying stars, which maintains the gas density always above the threshold level 
and which produces SF histories substantially different from those obtained with the 
standard IMF, for which the effect of the threshold is remarkable, in particualr in the 
last 4 Gyr of evolution. On the other hand, the case with /? = 2.35 shows a dependence 
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on the SF threshold even larger than the standard case, and this is due to the lower mass 
ejection rates from dying stars which stem from a steeper IGIMF (see Fig. \Q. 
In the middle panel, we show the abundance pattern predicted for the three values of /?, 
compared to those predicted in the standard case. The elements studied here are 0,Si 
and Fe, since the theoretical understanding of their production is quite clear, and their 
study allows us to neglect uncertainties related to their nucleosynthesis. In general, 
lower values of yS produce higher [Fe/H] values and higher [or/Fe] values (i.e. more 
a-enhanced elemental abundances) at a given metallicity. This is related to the fact that 
the fraction of massive stars increases with decreasing /3 and considering that massive 
stars are the main producers of a elements. Similar conclusions can be drawn by look- 
ing at the plots of the stellar metallicity distributions: the larger the value of yS, the larger 
the relative fraction of stars producing Fe, i.e. mostly type la SNe, i.e. stars in binary 
systems with initial mass ranging from O.SMq to SM©, hence larger the Fe abundances 
at any given epoch. This translates in SMDs peaking at higher [Fe/H] values for lower 
values of y6, assuming the same SF efficiency. 

As can be seen from the right panel of Fig. \T\ the model calculated with the IGIMF pro- 
viding the best results is the one with (3 = 2 and SF efficiency 0.3-0.5 Gyr"^ . The results 
obtained with this choice of yS are quite similar to the ones achieved with the standard 
IMF. This should not be a surprise since, as shown in Fig. [T] in the intermediate case 
with P = 2 the IGIMF is very similar to the standard IMF. The assumption of p = 2 
allows us to satisfactorily reproduce the set of observational constraints considered in 
this work. This is an important result, given the fact that the IGIMF is computed from 
first principles. On the basis of the results described in this section, it may be difficut 
to discriminate between the scenario with the standard IMF and the IGIMF with yS = 2. 
In the next section, the use of a diagnostic possibly useful to disentangle between the 
standard IMF and the IGIMF will be discussed. 



5. Discussion 

In this paper. We have modelled the physical properties of the S. N. within the IGIMF 
theory. In this scenario, the IGIMF can be calculated by combining the cluster mass 
function with the stellar IMF, which represents the mass function of stars bom within 
clusters and which can be described by a double-slope power law. An important feature 
of the IGIMF is that it depends on the star formation rate, which in turn evolves with 
time. 

The parameter /? regulating the cluster mass function may have an important im- 
pact on the predicted properties of the Solar Neighbourhood. In general alower value 
for y6 corresponds to a flatter IGIMF. In terms of chemical evolution, a flatter IGIMF 
translates into higher mass ejection rates from dying stars, hence globally a lower mass 
fraction incorporated into stellar remnants and higher gas mass densities. This implies 
that the evolution of all the models computed assuming (3=1 and p = 2 aie not sen- 
sitive to the star formation threshold and the star formation histories do not exhibit the 
"gasping" features typical of the standard model, which in turn is dominated by thresh- 
old effects at evolutionary times greater than 10 Gyr. Moreover, the lower the value of 
P, the higher the SN rate, the higher the metallicity and the larger the a-enhancement 
visible in the abundance pattern. The statistical test used to compare model results and 
the obervables indicates that the model which best reproduces the local observables is 
carachterized by /J = 2 as the index of the CMF The results of the best model are very 
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similar to those obtained with the standard case. A possible diagnostic which could 
help us disentangling between the two is represented by the present-day mass function 
(PDMF). The PDMF represents the mass function of living stars as observed in the so- 
lar neighbourhood. This quantity is an important diagnostic since it provides pieces of 
information complementary to the ones from the previously discussed observables. 

In the left panel of Fig.[3j we show the PDMF observed in the S. N. and predicted 
by means of our models. The PDMF computed with the standard IMF agrees with the 
observations in the range 0.4 Mq - 2 Mq. At very low stellar masses, the standard IMF 
seems too steep, whereas the distribution of stars with masses > SOM© is underesti- 
mated. Once again, this is due to the SF threshold, which has strong effects on the SF 
history of the solar neighbourhood at late times, inhibiting recent SF and hence causing 
the underabundance or absence of very massive stars. In contrast, the models calculated 
with the IGIMF provide all similarly a very good fit to the observed PDMF. 
The analysis of Fig. [3] seems to suggest that the SF threshold should not play a dominant 
role in the late evolution of the S. N. Within the IGIMF theory, the existence of a SF 
threshold may be an observational selection effect, naturally explained in this context as 
shown by Pflamm-Altenburg et al. (these proceedings). However, as shown by chem- 
ical evolution results, the SF threshold is fundamental in reproducing the metallicity 
gradients observed in the MW and in local galaxies, unless a variable star formation ef- 
ficiency through the disc is assumed (Colavitti et al. 2009). The study of the abundance 
gradients within the IGIMF theory may be of crucial help in sheding light on this issue 
and will be considered in future work. 

Another important issue concerns the time evolution of the IGIMF. In the right panel 
of Fig. [3l we show how the IGIMF varies as a function of time in the case of the three 
best models computed with different values of p. The best model (Ji = 2) shows very 
small variation of the IGIMF with cosmic time. Strong variations are predicted by as- 
suming j6 = 2.35, since in this case the star formation history is very much influenced 
by the effects of the SF threshold. It may be interesting to test the effects of the IGIMF 
in spiral. Milky Way-like galaxies in a cosmological context. Cosmological semiana- 
lytical models predic strong variations in the star formation histories of spiral galaxies 
(Calura & Menci 2009), which present a large number of spikes due to merging events 
and which should manifest into strong variations of the IGIMF with redshift. 
The universality of the IGIMF is another issue that deserves particular attention in the 
future. A chemical evolution study of elliptical galaxies within the IGIMF theory shows 
that the best value for j3 is 2.35, allowing to reproduce best the integrated a/Fe ratios 
observed in the local early-type galaxies. This value is in contrast with the best value 
suggested by the analysis of the S. N. features. A further study of the IGIMF in local 
dwarf irregular galaxies and dwarf spheroidals could certainly be of some help in this 
regard. 

Currently, the slope of the IGIMF in extreme SF conditions is another largely debated 
topic. Various indirect indications in external galaxies (Dabringhausen, these proceed- 
ings) and in our Galaxy (Stolte, these proceedings) seem to suggest that in strongly star 
forming systems, the slope of the stellar IMF should be flatter than the Salpeter one. 
Moreover, the assumption of a slightly top-heavy IMF in starbursts helps alleviating 
the discrepancy between cosmological models and observations regarding the a/Fe-cr 
relation observed in local ellipticals (Calura & Menci 2009). Addressing this subject 
within the IGIMF theory will be of primary importance in the nearest future. 
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Figure 3. Left panel: present-day mass function as predicted by means of our 
models compared to local observations. Right panel: time evolution of the IGIMF 
assuming various values for the parameter fi. 
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